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We initiate the development of a theory of the elasticity of nanoscale objects based upon new 
physical concepts which remain properly defined on the nanoscale. This theory provides a 
powerful way of understanding nanoscale elasticity in terms of local group contributions and 
gives insight into the breakdown of standard continuum relations. We also give two applications. 
In the first, we show how to use the theory to derive a new relation between the bending and 
stretching properties of nanomechanical resonators and to prove that it is much more accurate 
than the continuum-based relations currently employed in present experimental analyses. In 
the second, we use the new approach to link features of the underlining electronic structure to 
the elastic response of a silicon nanoresonator. 
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I INTRODUCTION 



The recent development of artificial free-standing 
structures of nanometer dimensions has led to great in- 
terest in their mechanical properties. A wealth of [experi- 
mental information is now available for nanowiresEHa and 
nanotubesJlEM and a. computational literature is devel- 
oping on the subject .ErtLij Many of these works make use 
of results from the continuum theory of elasticity to ana- 
lyze the behavior of nanometer structures. However, the 
applicability of continuum theories to nanoscale objects, 
where atomic-level inhomogeneities come to the fore, has 
yet to be explored in depth. 

Rigorous understanding of the elastic properties of 
nanoscale systems is crucial in understanding their me- 
chanical behavior and presents an intriguing theoretical 
challenge lying at the cross-over between the atomic level 
and the continuum. In the absence of an appropriate the- 
oretical description at this cross-over, critical questions 
remain to be answered including the extent to which con- 
tinuum theories can be pushed into the nanoregime, how 
to provide systematic corrections to continuum theory, 
what effects do different bonding arrangements have on 
elastic response, and what signatures in the electronic 
structure correlate with the mechanical properties of the 
overall structure? 

Recently, there have been a number of theoretical ex- 
plorations of the imnact of nanoscale structure on me- 
chanical properties.Ejlia These studies fall under two 
broad approaches, either the addition of-suxface and edge 
corrections to bulk continuum theoriescJ'EJ or the extrac- 
tion of overall mechanical response from atomic scale in- 
teractions. E3li3 The latter approach has the distinct ad- 
vantage of allowing first principles understanding of how 
different chemical groups and bonding arrangements con- 
tribute to overall elastic response, thus opening the po- 
tential for the rational design of nanostructures with spe- 
cific properties. 

In coarse graining from interatomic interactions to me- 



chanical response, some works rely upon the problem- 
atic decomposition of the total .system energy into a di- 
rect sum of atomic energies ,E3"Ej which is always arbi- 
trary and particularly inconvenient for connection with 
ab initio electronic structure calculations. The remain- 
ing works which attempt to build, up overall response 
from atomic level contributionsEZrllj fail to account prop- 
erly for the Poisson effect. Below we show that failure 
to account for this effect leads to surprisingly unphysical 
results. 

This manuscript presents the first theory for the analy- 
sis of overall mechanical response in terms of atomic- level 
observables which suffers from neither difficulty from the 
preceding paragraph. This analysis allows, for the first 
time, quantitative understanding of how continuum the- 
ory breaks down on the nanoscale, of how to make appro- 
priate corrections, and of how to predict the effects of dif- 
ferent bonding arrangements on overall elastic response. 
It is well known that the decomposition of overall elastic 
response into a sum of atomic level contributions is not 
unique. We show here, however, that with the additional 
constraint of dependence of moduli on local environment 
our definition of atomic level moduli becomes physically 
meaningful and essentially unique when coarse-grained 
over regions of extent comparable to the decay range of 
the force-constant matrix. 

For concreteness, in this work we focus on nanowires. 
However, we will also describe briefly how this work can 
be extended to any system with nanometer dimensions. 
The manuscript proceeds as follows. Section || briefly 
overviews the present state of the field. Next, as the tra- 
ditional concept of Young's modulus becomes ill-defined 
on the nanoscale, we begin by carefully defining contin- 
uum elastic constants appropriate for nanowires in Sec- 
tion [II . We then show how to decompose these constants 
exactly into atomic-level contributions based on true 
physical observables (rather than individual atomic en- 
ergies) using a straight-forward-application of Born and 
Huang's method of long wavesJ13 resulting in a decompo- 



sition similar in spirit to those in References [17] and [18] 
(Section IV). Section [v| demonstrates the surprising, rad- 
ical breakdown of this approach when applied to nanores- 
onators. Then, in Section fv|, we identify the source 
of the difficulty as the Poisson effect and present the 
first analysis of mechanical response truly applicable to 
nanoresonators. 

T he m anuscript then goes on to applications. Sec- 
tion VII verifies the physical meaningfulncss of our newly 



defined quantities by verifying that they predict response 
to modes of strain for which t hey w ere not directly con- 
structed. We then, in Section VIII, use our approach to 



generate a new, much more accurate, relationship be- 
tween experimentally accessible observables describing 
response to flexural and extensional strain in nanomc- 
chanical resonators. Finally, Section IX uses this theory 



to explore possible links between underlying electronic 
structure and local elastic response. 



II OVERVIEW 

As the introduction mentions, the literature pursues 
two broad categories of approach to the study of mechan- 
ical properties on the nanoscale, either surface and edge 
corrections to continuum theory or extraction of over- 
all response from the underlying atomic interactions. In 
the former category, Reference [13], through scaling ar- 
guments and numerical examples, notes that the Young's 
modulus for nanomechanical resonators scales as a bulk 
term plus surface and edge corrections. Although provid- 
ing insight and motivation, this work leaves completely 
open how one should understand these corrections from 
first principles. Reference [14] provides a more rigorous 
study based on separating nanoscale systems into con- 
tinuum surface and bulk regions. This latter approach 
allows prediction of changes in stiffness properties as one 
approaches nanometer length scales and has the appeal of 
generating physically motivated correction terms. How- 
ever, it relies on the separation of a nanomechanical res- 
onator into bulk and surface continua as an ansatz and 
therefore neither predicts when such a picture suffices to 
give an accurate description nor prescribes further cor- 
rections. 

References [15-18], on the other hand, start from the 
more general atomic level description and then try to un- 
veil physical properties from the underlying atomic de- 
scription. It is important to note that these works do not 
deal directly with nanoscale systems but rather focus on 
the effects of nanoscale inhomogeneities in bulk systems. 

References [15] and [16] concern the elastic proper- 
ties of grain boundaries. These works define atomic-level 
elastic moduli as the second derivative of the energy as- 
sociated with each atom with respect to strain and then 
go on to study the behavior of such moduli near grain 
boundaries. The difficulty with this approach is that it 
requires a breakdown into individual contributions from 



each atom of the total energy of any system. Such an 
atomic energy is neither observable nor uniquely defined 
and therefore cannot serve as an appropriate basis for 
theoretical understanding. 

Although References [17] and [18] work from valid 
physical observables, the components of the force- 
constant matrix, these works focus on bulk-like or meso- 
scopic scale systems and fail for nanoscale systems for 
the reasons which we describe in this work. Reference 
[17] investigates nonlocal elastic constants on the meso- 
scopic scale and links them to the underlying atomic in- 
teractions. It then proceeds to define an elastic constant 
for each atom and studies the behavior of these quanti- 
ties near surfaces and grain boundaries. Reference [18] 
defines a bond frequency from the force-constant matrix 
from which it deduces the possibility of bond rupture dur- 
ing crack nucleation. Neither of the above works prop- 
erly accounts for the Poisson effect, which we show in 
Sections |v|, and |v| to play a critical role in the elas- 
ticity of nanoscale systems. Moreover, straightforward 
generalization of these works to include this effect fail for 
the the same reasons as does the related approach which 
we describe in Section |y|. 



Ill NANOWIRE RIGIDITIES 



The prime difficulty in the application of continuum 
theory to objects of nanometer cross-section is the loss of 
the ratio of the inter-atomic spacing to the cross-sectional 
dimension as a small parameter. However, so long as the 
length of an object and the wavelength of the distortions 
considered both greatly exceed the inter-atomic spacing 
and the cross-sectional dimension, the object properly 
may be viewed as a one-dimensional continuum. Al- 
though we focus in this work on nanowires, the gener- 
alization of the discussion below to nanoscale systems of 
other dimensionality such as thin plates or nanoscopic 
objects is straightforward. 

Viewed as a linear continuum, the free energy per unit 
length / of a nanowire is 



/= {Eu 2 + FR-' 2 + Tt 2 )/2, 



(1) 



where u is the linear strain of extension, R is the ra- 
dius of curvature and r is the rate of twist of the tor- 
sion. The coupling constant E is the extensional rigidity, 
F is the flexural rigidity and T is the torsional rigid- 
ity. Unlike traditional bulk continuum concepts, the free- 
energy function Eq. (|l]) is observable in principle and 
thus provides an unambiguous operational definition of 
the rigidities. We avoid the use of traditional contin- 
uum concepts, such as the Young's modulus and the 
cross-sectional area, because such concepts are neither 
uniquely nor well-defined for nanoscale systems. 
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The rigidities in Eq. ([[]) are related to the phonon fre- 
quencies through 



LU TA 



= y/E/(\m)q, 
= v/F/(Am)g 2 , 



= y/T/(\mI r )q, 



(2) 
(3) 
(4) 



where to is the frequency for either the longitudinal, 
transverse or rotational acoustic modes, respectively, A 
is the linear atomic number density, q is the wave vector 
and m is the mass of a single atom. (This work focuses 
on single species systems for simplicity.) Finally, I r is 
defined unambiguously as the mean rotational moment 
I r = {l/N c ) J2 a ( x a + Va)' where the sum ranges over all 
atoms in the cell, N c is the number of atoms per unit 
cell, the wire is assumed to run along the z-axis and the 
origin lies on the center line of the wire. 

Finally we note that although the rigidities in Eq. (Q) 
are well-defined, certain traditional continuum relations 
between them do not hold. Specifically, we will show be- 
low that the traditional continuum relationship between 
the extensional rigidity E and the flexural rigidity F fails 
on the nanoscale, similar to relations whicrU^ws-been re- 
cently used in the analysis of experimentsJiiBE 



IV METHOD OF LONG WAVES 

One reason for breakdown of traditional continuum re- 
lations on the nanoscale is that the continuum perspec- 
tive course grains away important fluctuations which oc- 
cur over distances on the order of the inter-atomic spac- 
ing. To overcome this shortcoming, we propose to coarse 
grain only on distances over which the underlying inter- 
atomic interactions vary, the decay length of the force- 
constant matrix. The straightforward approach to gen- 
erate such a theory is the '.iaethod of long waves" devel- 
oped by Born and Huang fB which is somewhat similar 
to the approaches whiclt-have been used previously to 
defects in bulk systems £3' Ea This section applies the 
method of long waves to nanoresonators. The next sec- 
tion shows how, surprisingly, this and related approaches 
fail in the study of systems with free surfaces and there- 
fore, in general, cannot be used to describe nanoscale 
systems. In Section VI, we describe how to go beyond 
the straightforward application of the "method of long 
waves" in order to achieve a meaningful description. 

We focus initially on longitudinal waves and, as noted 
in Section |lj], choose our origin to lie on the center 
line of the wire, which we let run along the z-axis. For 
nanowires, the presence of surfaces breaks periodicity in 
the transverse directions resulting in a one-dimensional 
crystal with an extremely large unit cell of length L c = 
N c /X, where N c and A are as in Section III. In all expres- 
sions below, boldfaced quantities are 3-ZV c -dimensional 
and arrowed vector quantities are three-dimensional. Fi- 
nally, sums with Greek indices range over atoms in the 
unit cell. 



To relate the rigidities to the dynamical matrix, we 
begin similarly to Born and Huang and choose to factor 
the Bloch phases (e lqz ) out of the representation of the 
phonon polarization vector u, incorporating them into 
the definition of the dynamical matrix D, so that the 
acoustic phonon polarization vectors are periodic across 
the cell boundaries. This ensures a uniform description 
of the distribution of elastic energy along the axis of the 
wire. To generate a scalar equation for the phonon fre- 
quency, Born and Huang project the secular equation for 
the dynamical matrix, 



Du 



-mu> 2 u, 



(5) 



against the zeroth-order polarization vector . Here, 
however, to more symmetrically represent the distribu- 
tion of elastic energy, we project against the full polar- 
ization vector u. Equating the frequency lu in Eq. (|5|) 
with the longitudinal frequency in Eq. (H) gives 



E _ [utDupl 
A utu 



£ — 

s,t=0 



u tW D P 



(6) 



where we have expanded the numerator of the Rayleigh 
quotient to second-order in powers of (iq) and where the 
3x3 sub-block of D^', which couples atoms a and /3, is 



I DN ]^ = ^ E (z-(R + r - f a )) n . 



(7) 



Here, R is a lattice vector along the z — axis, <5> a p(R) 
is the 3x3 sub-block of the force-constant matrix which 
couples atoms a and (3 located at positions f a and R+fp, 
respectively, and 



d 2 U 



dr a d(R + fp) 



Finally, substituting Eq. ((?]) and [u'°l] Q = z into 
Eq. (^), allows us to express E as a sum over atoms (a) 
in the unit cell and all atoms (/?, R) in the system, 



E 



L 

a 

PR 



AZ a p-<f> a p(R)-AZ af3 /2- 



(8) 



*a/j(i?) -uf 



+AZ al3 - <t> al3 (R) ■ uf - fiM • $ a/3 (R) -AZ af3 }. (9) 

Here, AZ a p = zz ■ (f a — fp — R), !?„ is the first-order 
polarization vector and we refer to the e a as the "atomic 
moduli" . 

The atomic moduli as currently defined in Eq. ([)]) pro- 
vide a useful microscopic analysis of elastic response in 
bulk systems which is similar in spirit to the decomposi- 
tions rH6p4 previously in the study of bulk material sys- 
temsB'El To see that Eq. (O) indeed decomposes the 
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overall elastic response of bulk systems into atomic con- 
tributions coarse-grained over distances on the order of 
the decay-length of the force-constant matrix, we note 
first that for infinite bulk systems, elastic waves are pla- 
nar. This implies that the first-order polarization vec- 
tor («„') is uniform from primitive cell to primitive cell 
and thus depends only on the local environment of each 
atom. Next, we note that although the strain terms 
(AZ a p) scale linearly with distance between atoms, the 
terms which contribute to the final result are bounded in 
range by the inter-atomic interactions ($(i?)). Thus, the 
atomic moduli depend only on the local atomic environ- 
ment over distances which the decay of the force-constant 
matrix determines. 



structure. Finally, we employ periodic boundary condi- 
tions along the ^-direction. 

Figure |2|a shows that the atomic moduli e a predicted 
for nanowires using the straightforward approach of 
Eq. (||) are unphysical in that they depend upon the 
macroscopic dimensions of the system and not simply 
on the local environment of each atom. In particular, 
the atomic moduli on the surface grow linearly with the 
diameter of the wire and the moduli in the center of the 
wire fail to approach the expected bulk limit, ej, = Yb/ pb, 
where Yf, is the Young's modulus in bulk and pb is the 
number density in bulk. (Note that these two effects are 
interelated, as the moduli must sum to give the macro- 
scopic value in the bulk limit.) 




FIG. 1: Atomic structure of silicon nanowires with an ap- 
proximate cross-sectional diameter of 1.5 nm. The wires 
are viewed along the longitudinal axis: c(2 x 2) structure 
(left), 2x1 structure (right). 



FAILURE OF STRAIGHTFORWARD 
APPROACH IN NANOWIRES 



We now demonstrate through direct calculations that 
the approach outlined in the previous section gives un- 
physical results when applied to nanoresonators. Specif- 
ically, we study the behavior of [100]-oriented-Jaanores- 
onators of silicon, which recent ab initio studiesE3 predict 
to undergo a size-dependent structural phase transition 
between the two structures in Figure [l] at a cross-section 
of ^3 nm. (The interested reader may refer to Refer- 
ence [21] for explicit details of the microscopic structure 
of these wires.) Initially, we— work with the Stillinger- 
Weber inter-atomic potential^! which suffices for the ex- 
ploration of general nanoelastic phenomena and which 
allows study of cells with many thousands of atoms. 



Later in the manuscript (Section [X) we use the Sawada 
tight-binding mode£3 with modifications proposed by 
Kohyamao to explore the correlation between our local 
approach and the underlying electronic structure. For all 
calculations below, we fully relax the atomic coordinates, 
the periodicity of the wire and, need be, the electronic 
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FIG. 2: Predictions of atomic moduli e a for c(2 x 2) 
nanowires of varying diameter: (a) straightforward the- 
ory (Eq. (§)) and (b) new theory (Eq. (|l])). The insets 
denote the approximate diameters of the wires. The value 
of the atomic modulus (Eq. (||) or Eq. ([ll])) is along the 
ordinate and radial distance of the atom from the center 
line is along the abscissa. 

Although one has some freedom in choosing the terms 
used in the perturbation expansion Eqs. (^|) and (^), for 
example to project Eq. (||) against u^l instead of u, all 
such expansions will lead to similar linear scaling along 
the surface of the wire and approach an incorrect value 
at the center of the wire. Thus, straightforward applica- 
tion of the method of long waves fails to result in a local, 
and hence physically meaningful, description of elastic re- 
sponse in nanoscale systems, as will straightforward vari- 
ations thereon such as those in References [17] and [18]. 



VI METHOD OF LONG WAVES IN 
NANOWIRES 

To cure the difficulties uncovered in the previous sec- 
tion, we proceed by first identifying the cause of the 
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pathological behavior and then exploiting the freedom 
in Eq. (|9|) to remove this pathology. 

The failure of the straightforward approach arises from 
the fact that elastic waves in nanoresonators, or any sys- 
tem with free surfaces, are not strictly planar. In partic- 
ular, the Poisson effect, which the first order polarization 
vector uM contains, causes each atom to displace by an 
amount in direct proportion to its distance from the cen- 
ter line of the system. Eq. @ then leads directly to linear 
scaling of the atomic moduli at the surfaces of the wire. 

Defining an atomic elastic reponse dependent soley on 
the local enivornment requires separation of extensive 
elastic effects from intensive nanoscopic effects. To seper- 
ate the extenive motion in the first order polarization 
vector from that of its intensive motion, we define the 
atomic displacements (u£f) as the intensive nanoscopic 
motions, 



djk<Jj, is a diagonal 3x3 matrix with el- 



it:! 



{-a x xx ■ -a y yy-)T a , 



(10) 



where a x>y are the Poisson ratios. There are three logical 
choices for the Poisson ratios: some sort of local atomic 
definition, an overall average for the wire, or the bulk 
values. We choose to use bulk Poisson ratios for a num- 
ber of reasons. First, a locally varying definition makes it 
impossible to exploit the continuous rotational and trans- 
lational symmetries in the dynamical matrix, which we 
find necessary to employ below in constructing atomic 
moduli with local behavior. Second, only by employing 
the bulk (rather than average) Poisson ratios do we find 
a definition which approaches the appropriate bulk value 
in the centers of wires of finite width. Finally, we note 
that we always have the freedom of working with bulk 
Poisson ratios because any motion along the surface in 
addition to that resulting from bulk Poisson effect will 
not scale extensively with the diameter of the wire and 
can therefore be incorporated into the intensive atomic 
displacements u£. 



rfrl 



After making the decomposition in Eq. (|1C|), we next 
employ the continuous rotational and translational sym- 
metries of the dynamical matrix to eliminate all extensive 
dependencies in Eq. (|J). Appendix ^] outlines the proce- 
dure for doing this, which then transforms Eq. (^) into 

e a = ^2{~ Af Q/3 ■ V a p(R) ■ Af Q/3 /2 + u"J ■ ® a0 (R) ■ 

PR 

+Ar al3 ■ $ a/3 (£) • - u r J ■ $ a p(R) ■ Af Q/3 }, (11) 

where Ar a p represents the total strain between atoms a 
and 0, 

Af a/3 = (-cJz xx - a y yy + zz) • (f a -fp-R), 
and ^ a p(R) renormalizes as 

¥ a p(ft) = 2$ Q/3 (i?) - (Tr^ af3 (R) ■ S) E~\ 



where 

ements o~ x , o~ y and a z = —1, respectively and 5jk is the 
Kronecker delta. 

This new construction ensures that the modulus of 
each atom depends only upon its local atomic environ- 
ment because u£ no longer includes extensive motions 
and, although Ar Q/ 3 still depends on relative atomic dis- 
tances, the renormalized \P decays as $ does. Thus, it is 
now the range of the force-constant matrix which controls 
the size of the neighborhood upon which each atomic 
modulus can depend. Therefore, the moduli of atoms 
in the interior now must correspond to the expected bulk 
value, the moduli of the atoms on the surface now cannot 
depend upon the extent of the system, and the resulting 
description is physically meaningful. Figure illustrates 
the success of Eq. ( |Tl| ) . 

The fact that decomposition of elastic response into 
atomic level contributions is not unique raises questions 
as to the physical meaning of such a decomposition. The 
new decomposition Eq. ( pd| ) is the first which remains 
dependent only upon local environment for systems with 
free surfaces. Any other definition which respects locality 
can only redistribute portions of each atom's modulus 
among other atoms within a region of extent compara- 
ble to the decay of the force-constant matrix. Any sum 
over such a region of the moduli will always be nearly the 
same. Therefore, coarse grained over such regions, prop- 
erly localized a tomi c moduli become physically mean- 
ingful. Section VII and Appendix demonstrate this 



explicitly by comparing the predictions for flexion from 
either properly localized moduli or straightforwardly de- 
fined moduli, respectively. To further demonstrate that 
alternate local definitions are equivalent in this coarse- 
grained sense, we have explored alternate local construc- 
tions. In particular, while our present construction takes 
care to employ the continuous symmetries of the dynami- 
cal matrix in such a way so as to respect symmetry among 
the x, y, z Cartesian coordinates, we have also repeated 
the construction while treating the x, y coordinates sym- 
metrically but not the z coordinate and have found nearly 
identical results for all of the applications below. 

We close this section with a brief description of how 
the above approach extends to any system of nanome- 
ter dimensions. To be considered small in this context, 
a dimension must be much smaller than the wavelength 
of the distortions considered. This manuscript focuses on 
nanowires, systems with two small dimensions. For a sys- 
tem with one small dimension, for instance a plate with 
nanometer thickness, the above approach develops in the 
same way with the one minor change that the definition 
of u£ (Eq. ( |io| ) involves the Poisson effect only in the one 
small dimension. For an object with three small dimen- 
sions, for instance adiabatic loading of a nano-object, the 
nature of the Poisson effect depends upon the mode of 
loading, and the system should treated as either of the 
above cases accordingly. Finally, objects of no small di- 
mension, and hence no Poisson effect, can be described 
straightforwardly as bulk-like using Eq. (O). 
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VII TRANSFERABILITY 

To establish that our new atomic moduli are not mere 
convenient mathematical constructions but are physi- 
cally meaningful, we now consider their transferability 
to phenomena not considered in their original construc- 
tion. In particular, we consider flexion, where the elas- 
tic distortion is no longer homogeneous throughout the 
cross-section. 

If our atomic moduli are indeed a measure of the local 
elastic response, then under flexion the free energy per 
unit length will take the form 



(a) 



/ = 



(l/L c )£>al£/2, 



(12) 



where u a is a measure of the longitudinal strain which 
atom a experiences. For this form to be sensible, the 
diameter D of the wire must not become comparable to 
the range of the force-constant matrix so that the atoms 
which contribute to each e a all experience similar strains 
u a . 

Within continuum theory, uniform flexion with radius 
of curvature R corresponds to a longitudinal strain which 
varies linearly across the wire, u = x/R. (This holds to 
better than to two parts in 10 3 for all wires in our study.) 
This would then predict a flexural rigidity of 



Fat — j ^ y Ga^Cct 



(13) 



Figure |^a shows the fractional error (SF at ) 

SF at ee (F at - F)/F, (14) 



in predicting the flexural rigidity from Eq. (13), where 
F is determined directly through numerical calculations. 
The figure shows that these errors are indeed quite small. 
Note that use of the straightforward definition in Eq. (|]) 
with its unusually scaling surface moduli leads to invalid 
predictions for flexion. (Appendix [b| shows this directly 
through scaling arguments.) As a result, to have predic- 
tive power, definitions of atomic moduli must properly 
account for the Poisson effect. 

To demonstrate that our new approach has greater pre- 
dictive power than traditional continuum approaches, the 
figure also shows the fractional error (SF tc ) in predicting 
the flexural rigidity when using the traditional continuum 
relation 



Ftc = E(I/A), 
SF tc = {F tc ~F)/F 



(15) 
(16) 



where I /A defines the mean bending moment, which we 
define unambiguously as 



1/ A = (1/A C )]T; 



(17) 
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FIG. 3: Fractional error as a function of inverse diameter 
1/D in predicting (a) flexural rigidity from extensional 
properties (SF at , Eq. (|l|) and SF tc , Eq.Q) and (b) ex- 
tensional rigidity from flexural properties (SE at , Eq. j2fj| ) 
and SE tc , Eq. (|2l|)). Note, for convenience, the abscissa 
is labeled by D and not 1/D. 



The linear behavior of the fractional error SF tc as a 
function of 1 / D indicates that continuum theory does not 
properly account for surface effects, a result of the fact 
that flexion places a larger emphasis on the surface than 
does extension. The dramatic improvement from the use 
of Eq. ( |l3| ) arises because the atomic moduli place proper 
emphasis on the surface and on the interior, as they prop- 
erly treat each atomic environment locally. The atomic 
moduli therefore properly account for elastic fluctuations 
along the cross-section of the wire which are on scales too 
small for traditional continuum theories to capture. 

Finally, we note that the only appreciable error within 
the new framework occurs for the smallest wires (Dps 
1.5nm). At this point, the cross-sectional dimension be- 
comes comparable to the range of the force-constant ma- 
trix, and Eq. ( p"2| ) represents an improper use of the phys- 
ical concept of atomic moduli. As described above (sec- 
ond from last paragraph in Section VI), only sums of 



atomic moduli over regions of extent comparable to the 
range of the force-constant matrix carry physical mean- 
ing. Any sum sensitive to variations over shorter scales, 
as is Eq. (12) when limited to wires narrower than the 
range of the force-constant matrix, cannot be depended 
upon to lead to meaningful results. This underscores 
the fact that properly construed atomic moduli are not 
truly atomic-level quantities but a concept coarse grained 
over the range of the force constant matrix. As we have 
seen, however, this coarse-graining is on scales signifi- 
cantly smaller than those captured by traditional contin- 
uum theory. 
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VIII EXTENSION FROM FLEXION 

We now apply the above denned atomic moduli to de- 
rive a new relation for the extensional modulus in terms 
of the flexural modulus and other experimentally accessi- 
ble obserjvahies, a relation often needed in experimental 
analyses.BB&El We then show that the new relation is 
much more accurate than the standard continuum-theory 
based relation currently employed in experimental anal- 
yses. Finally, we employ our concept of atomic moduli 
to provide quantitative insight into the improvement of 
our new relation over the traditional continuum relation. 

The basis for the following analysis is the fact that 
Eq. (13) gives a very good estimate of the true flexural 



modulus, as Figure 
for E, Eqs. (0) anc 



ia confirms. Using the exact relation 
(|ll|), and Eq. (p^|), one derives the 
following leading-order prediction E at for the extensional 
modulus, 



1 / F A 
= 2\JjA +Yb J b 

with a predicted error SE^ of 



(18) 



8E$> = 



EL, 



{B at [(e s ) x 2 - e b ] + {(e s ) x 2 - (e s )]}(19) 



Here, N s is the number of "surface" atoms, defined as 
those for which e a differs significantly from the limiting 
bulk value ej,, (e s ) x 2 is the inertia weighted average sur- 
face moduli 



(e s ) 3 



_ Es e s x s 



2 

S "~8 



with the sums (E s ) ranging over "surface" atoms, (e s 
is the average surface moduli 



— V, 



and 



B„ 



(1/2) 



iV a;2 

J_V x 2 



where ^ c implies sums over all atoms in the unit cell. 
This result holds for any division of the atoms into "sur- 
face" and "bulk" to the extent that each "bulk" atom 
has atomic modulus e b - 

Figure gb shows the relative error SE at , 



SE at = {E at - E)/E, 



(20) 



between the extensional modulus E determined directly 
from numerical calculation and as determined from our 
new relation, Eq. (|lg|). (Note that SE^ — SE at exactly 
when the moduli prediction F at = F holds.) As the 



relevant point of comparison, the figure also shows the 
relative error 5E tc , 



6E tc = (E tc - E)/E, 

associated with the traditional continuum result 

,.4 



E, 



F- 



(21) 



(22) 



Hence, Eq. (|18|) is much more accurate than the standard 
continuum result Eq. (|22|). 

To understand the improvement of the new relation, 
Eqs. (§) and @ may also be combined to yield a predic- 
tion for the fractional error in the traditional continuum 
analysis, 



5E\ 



(p) 



EL C 



{B tc [(e s ) x 2 - e b ] + [{e s ) x 2 - (e s )]} , (23) 



which takes precisely the same form as Eq. ( p^ ) except 
for the change in the prefactor in the first term from B at 
to B tc = 2B at + 1 . For continuous wires of homogeneous 
circular or regular polygonal cross-section, we have ex- 
actly B at = 0. Thus, generally we expect B at to be close 
to zero and B tc to be close to unity. We now note that 
the term in the first set of square brackets ([(e s ) x 2 — e ]) 
in both Eq. ( fl9| ) and Eq. (23) is an average difference 
between surface and bulk atoms and, therefore, is gen- 
erally much larger than the term in the second set of 
square brackets [{e s ) x 2 — (e s )], which is the difference be- 
tween two differently weighted averages over the surface 
atoms. Thus the larger term nearly vanishes in our new 
relation, Eq. (|19|) , but not in the traditional continuum 
relation, Eq. (|2^). From this analysis, we see that the 
reason why the traditional continuum relation has larger 
errors is that it does not properly differentiate between 
the local surface and bulk environments. 

The atomic moduli also lead to a quick, intuitive ar- 
gument to understand the improvement of Eq. ( p"8| ) over 
Eq. ( |22| ) . From our results we know that the surface mod- 
uli can be quite different than those of the bulk. It is also 
known that flexion places larger emphasis on the surface 
than does extension. If the average surface modulus is 
less/more than that of the bulk then the first term in 
Eq. ( |l8| ) (F ■ A/ 1) will underestimate/overestimate the 
extensional rigidity, while the second term will overes- 
timate/underestimate it. Therefore, errors will tend to 
cancel in the average of the two. 



IX CORRELATION BETWEEN LOCAL 
ATOMIC MODULI AND ELECTRONIC 
STRUCTURE 

We now explore the local atomic moduli as a link be- 
tween local elastic properties and the underlying elec- 
tronic structure. To do this, we have calculated the 

sth the c(2 x 2) and (2 x 1) 
1 tight-binding model with 



atomic moduli, Eq. (11), for 
structures, using the Sawada 
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modifications proposed by Kohyama.Eil We have stud- 
ied various wires, all of which give similar results. For 
brevity, we here only report on wires with cross-sectional 
diameter rj 2.4nm. The supercell of our calculation is 
four bulk cubic lattice constants long in the periodic di- 
rection, and hence sampling the Brillouin zone at the T 
point in the electronic structure calculations is more than 
sufficient. 

The left panels of Figures |^ and || show the resulting 
atomic moduli, with values color coded so that yellow cor- 
responds to the bulk value (17.1 eV/atom). The figure 
shows that deviations from this value concentrate near 
the surface in patterns characteristic of the structure of 
the wire. Moduli near the surface fluctuate widely, rang- 
ing from 6-27 eV/atom for the c(2 x 2) structure and from 
4-26 eV/atom for the 2x1 structure. 



KM I I I* 

~ V' V V' V V v. 

C I. I I 11. I 
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FIG. 4: c(2 x 2) wire, with cross-section rj 2.4nm. Re- 
sults are from the tight-binding model described in the 
text. Left: atomic moduli color coded from small to 
large: blue-green- yellow-red. Right: valence electronic 
density projected (integrated along the wire axis) onto 
the cross-section color coded from low to high: blue- 
green-yellow-red. White dots in the right figure indicate 
the location of atom cores. Note the correlation of large 
atomic moduli with large charge density along the sur- 
face. 




FIG. 5: 2 x 1 wire, with cross-section 
conventions as Figure 0. 



2.4nm: same 



To allow comparison with the underlying electronic 
structure, the right panels in the figures display the va- 
lence charge densities from the tight-binding calculation 
projected (integrated along the wire axis) onto the cross- 
section of the wire. (To compute the electron density 
from the tight-binding coefficients, we employed orbitals 



from a density functional calculation of the silicon atom.) 
The figures display the electron densities using a color 
map similar to that employed for the atomic moduli. In- 
triguingly, there is an apparent correlation between the 
values of the atomic moduli and the underlying electron 
density. In particular, large/small atomic moduli corre- 
late with regions of large/small electron density in Fig- 
ure [|, indicating that charge distribution along the sur- 
face of these wires greatly affects the local elastic proper- 
ties and thereby the overall elastic response, particularly 
to flexion which emphasizes surface effects. Figure ^| ex- 
hibits a similar correlation, but not as pronounced. Note 
that, in this figure, coincidence of red atomic moduli just 
under the surfaces of the first layer of atoms correlates 
with red charge densities in the same location. 

Unlike in the classical potential case, where properly 
defined surface moduli are systematically lower or equal 
to the bulk value due to decrease in the number of bonds 
(Figure |^b), we find that in a quantum model, surface 
moduli may even greatly exceed the bulk value due to 
changes in local charge density which can enhance the 
mechanical strength of bonds (Figures | and |.) This 
contrast underscores both the importance of considering 
contributions of the electronic structure to mechanical 
response and the need for a definition of atomic moduli 
which can be computed from physical obscrvablcs obtain- 
able from electronic structure calculations. 

The ultimate use of atomic level moduli is to under- 
stand mechanical response. Table [| compares the errors 
from both continuum theory and the use of atomic-level 
moduli in predicting mechanical response for the two 
wires under consideration in this section. The quanti- 
ties compared exactly parallel those of the previous sec- 
tion. The first two columns of the table consider predic- 
tion of flexural response from continuum theory and our 
atomic moduli, respectively, and the second two columns 
consider prediction of extensional response from flexu- 
ral response using either the traditional continuum rela- 
tion or our new relation, respectively. The table shows 
that the new relations, Eqs. ([H]) and ( |l8| ) are again 
very accurate. The table also shows that these predic- 
tions are superior to the corresponding continuum re- 
sults, Eqs. fll5|) and (22), respectively. 



Structure 


5F tc 


SFat 


5E tc 


SE at 


2 x 1 


13.1% 


2.86% 


-11.6% 


1.78% 


c(2 x 2) 


4.91% 


1.79% 


-4.68% 


2.02% 



TABLE I: Comparison of errors between the tradi- 
tional continuum theory and the atomic moduli de- 
scription when predicting flexural response through 
Eqs (H|) and dH) (fi rst and second columns, respectively) 
and when inferring extensional response from flexural 
response through Eqs. (^) and (fL8| ) (third and fourth 
columns, respectively). 

Interestingly, the continuum predictions for the c(2 x 2) 
wire are fairly reliable. The atomic moduli provide an av- 
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enue for understanding this as well. Fluctuations in the 
moduli in the c(2 x 2) wires are localized and hence av- 
erage out over regions of extent comparable to the decay 
of the force-constant matrix. Moreover, in this partic- 
ular case they tend to average to values close to that 
expected of the bulk. Without meaningful fluctuations 
on the length scales of the decay of the force-constant 
matrix, we expect continuum theory to perform well for 
this wire. In contrast, the 2x1 wire exhibits much more 
systematic variations in the moduli. The outermost sur- 
face atoms have a consistent and significantly reduced 
modulus, and there is also a clear significant and system- 
atic variation in the moduli throughout the cross-section 
of the wire. Because this second wire does exhibit mean- 
ingful fluctuation over distances comparable to the decay 
range of the force constant matrix, we expect traditional 
continuum relations to give particularly poor results, un- 
derscoring the importance of the local atomic-level mod- 
uli description. 



X CONCLUSIONS 

This manuscript presents the first definition of atomic- 
level elastic moduli for nanoscale systems which are de- 
fined in terms of physical observables, correctly sum to 
give the exact overall elastic response and depend only 
on the local environment of each atom. Although these 
moduli are not necessarily uniquely defined, their sum 
over regions of extent comparable to the range of the 
force constant matrix is physically meaningful and may 
be used to make accurate predictions of mechanical re- 
sponse. The moduli resulting from our formulation trans- 
fer to different modes of strain and correctly account for 
elastic fluctuations on the nanoscale. They also lead to 
a quantitative understanding of when traditional contin- 
uum relations breakdown and how to properly correct 
them properly. Specifically, we demonstrated a more ac- 
curate method for relating extensional and flexural prop- 
erties. These moduli provide a clear and natural method 
for distinguishing mechanically between "surface" atoms 
and "bulk" atoms and give insight into the correlation 
between the local mechanical response and the underly- 
ing electronic structure. Finally, these moduli allow the 
identification of which atomic arrangements lead to more 
pliant or stiffer response opening the possibility of their 
use as a tool to aid in the rational design of nanostruc- 
tures with specific mechanical properties. 
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A USE OF CONTINUOUS SYMMETRIES TO 
PRODUCE LOCAL MODULI 

This appendix outlines the use of rotational and trans- 
lation symmetries of the force constant matrix to refor- 
mulate the ill-defined atomic moduli in Eq. @ into the 
well-defined form in Eq. (p~T|) . 

From continuous translational symmetry, all force- 
constant matrices obey 

Y/b<*AR)-te = o, (i) 

S3,R 

for any vector cp — c that is the constant vector for all 
atoms in the unit cell. Moreover, the SO(3) rotational 
symmetries imply 

$ « g „3 s {R)rpt = $ ag< i3t (R)rf3 s , (2) 
/3,-R P-.R 

where the g, s, t correspond to one of the x, y, z Cartesian 
coordinates, $ a q,f3s(R) is the q, s component of & ai p(R) 
and rp s corresponds to the s component of the position 
vector rp. 

The ill-defined atomic moduli, Eqs. (|^)-(|9|), contain di- 
vergent terms which are in one of the following two forms: 

Y ^sr a s^oi S ,pt{R)u r p t , (3) 

a, PR 

or 

^ <y s r as ^asM{ R ) r 0t^t- (4) 

a,(SR 

All of these terms scale linearly along the surface of the 
wire and give rise to the linear scaling of the surface mod- 
uli with system size evident in Figure ^a. 
From Eq. (|l|), one can set Eq. (||) equal to 

Vs(r as ~ rp s )$ as ^ t (R)up t . (5) 

a,pR 

Using both Eqs. (||) and (||), one can set Eq. equal to 

- Y a * I ^ as ~~ r 0s)®as,f3t(R)(r a t ~ rpt) 
a, f3R 

-j(r as - rp s )§ at ^ t (R)(r as - rp s ) 

-\{r a t ~ rf3t)^ as ,ps(R){r a t - r 0t ) } &t- (6) 

The above terms now only depend on relative distances 
over a range controlled by the decay of the force-constant 







matrix and therefore no longer scale linearly with the 
size of the system. Combining the above transforma- 
tions with the separation of the extensive motion from 
the intensive motion in the first-order polarization vector, 
Eq. (RG), then results in the well defined form Eq. (O). 



B FAILURE OF STRAIGHTFORWARD 
APPROCHES TO PREDICT FLEXION 

We now demonstrate that although the moduli defined 
in Eq. (^|) indeed sum to give the correct overall exten- 
sional response, their unphysical scaling properties lead 
to invalid physical predictions when used in other con- 
texts and, therefore, that they are ill-defined. In particu- 
lar, we use simple scaling arguments to prove that these 
ill-defined moduli give a nonnegligible error in predicting 
the flexural rigidity in the continuum limit. 

Consider a circular wire whose radius R is sufficiently 
large such that continuum theory applies. From Fig- 
ure ||a we see that the surface moduli defined by Eq. (Q) 
scale linearly with R. Because of this, and in order to 
arrive at an analytic result for the error in predicting the 
flexural rigidity, we can assume, to a good approxima- 
tion, that all surface moduli e s are proportional to their 
radial distance r from the center of the wire, 

e s = e s r. 

Here, e s is the same for all surface atoms. Next, we 
define Ri as the inner radius such that all atoms at po- 
sition r < Ri are in the bulk and all atoms with po- 
sitions r > Ri are on the surface. Finally, we define 
AR = R-R t . 

Because of the facts that the surface moduli scale lin- 
early with the system and that the sum of all atomic 
moduli must equal the extensional rigidity we conclude 
that the moduli in the bulk region of the wire, derived 
from Eq. (||), e' b cannot equal the average value of the 
bulk material e&. (Figure ga also evidences this behav- 
ior.) The two above facts imply that the following equal- 
ity must hold, 

/•Ri r R 

e r dr, 

Jo 

or to leading order in AR/Ri, 

e s AR(l + AR/R^ + e' b /2 = e 6 (l/2 + AR/Ri). (1) 

In the continuum limit, the traditional continuum re- 
lation Eq. ([l5]) holds, and therefore the fractional error in 
predicting the flexural rigidity from the ill-defined atomic 
moduli Eq. (||) is equal to 



or to leading order in AR/Ri, 

e s AR(l + 2AR/Ri) + e' 6 /4 - e' b AR/R t 







/ e s r dr 4 


/ e' b rdr = 


'Ri 


Jo Jo 



SR 



ill — at — 



JiT de (fit, dr x2e * r + l R ' dr x2e 'b r ) 



SFill-at 



e b /A 



-(2) 



J d0 J Q dr x 2 ebr 



Using Eq. (|l|) to solve for e s Ai?, Eq. (||) becomes 
8F M - at = l-^+ 2(e6 ~ e ^ Afl/fl. 

Therefore, in the continuum limit, the ill-defined mod- 
uli do not approach the correct result and thus give a 
prediction which is even worse than that of traditional 
continuum theory. 
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